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CN \ Abstract 

CN ■ We study integrable and conformal boundary conditions for s£{2) Zjt parafermions 

. on a cylinder. These conformal field theories are realized as the continuum scaling 

, limit of critical A-D-E lattice models with negative spectral parameter. The confor- 

', mal boundary conditions labelled by {a,m) £ {G,1j2k) are identified with associated 

I integrable lattice boundary conditions labelled by (r, a) G {Ag-2,G) where g is the 

Coxeter number of the A-D-E graph G. We obtain analytically the boundary free 
energies, present general expressions for the parafermion cylinder partition functions 
and confirm these results by numerical calculations. 
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1 Introduction 

There has been much recent progress ||, H, |], |^ on understanding integrable bound- 
aries in statistical mechanics, conformal boundary conditions in rational conformal field 
theories and the intimate relations between them. Indeed it appears that, for certain classes 
of theories, all of the conformal boundary conditions can be realized as the continuum scaling 
limit of integrable boundary conditions for the associated integrable lattice models. For sim- 
ple rational conformal field theories, such as si{2) minimal theories, a complete classification 
has been given [|I], ^, |^ of the conformal boundary conditions. These are labelled by nodes 
(r, a) of a tensor product graph A ^ G where the pair of graphs {A, G), with G of A-D-E 
type, coincide precisely with the pairs in the A-D-E classification of Cappelli, Itzykson and 
Zuber [§]. Moreover, for the unitary minimal models the physical content of the boundary 
conditions has been ascertained 0] by studying the related integrable boundary conditions 
of the associated A-D-E lattice models p with positive spectral parameter u > 0. 

It is highly desirable to carry out similar studies on boundary conditions for other classes 
of conformal field theories and their associated lattice models such as the super conformal and 
higher fusion models, higher rank models and so on. In this paper we study the boundary 
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conditions of TL^ parafermionic models [IT], jT^], [13], |T1|, |T3|, jTB], |T^, [TB|, |TP[. This class 



of rational conformal field theories admits an A-D-E classification |TT|, similar to the 
classification of Cappelli, Itzykson and Zuber and is associated to the same A-D-E lattice 
models but in the regime with negative spectral parameter u < 0. The first two members of 
the A series are the Ising model (A3) and the hard hexagon model (^4) which has the 
same central charge as the 3-state Potts model but admits the W3 algebra as the extended 
chiral algebra and has 6 rather than 8 conformal boundary conditions. 

The layout of the paper is as follows. In Section we review the conformal properties 
of si{2) Zfc parafermions, their central charges, conformal weights, characters, Kac tables 
and modular matrices. We also discuss the modular invariant partition functions on the 
torus and the cylinder partition functions. In Section ^ we present the lattice realizations of 
these theories. We define the A-D-E lattice models and their integrable boundary weights, 
present their finite size corrections and compute their bulk and boundary free energies. In 
Section ^ we identify, through numerical analysis, the integrable boundary conditions from 
the lattice model with the conformal boundary conditions of parafermions. In Section ^, 
we discuss our numerical techniques and present some typical numerical data on which our 
results are based. We conclude in Section |] with a short discussion. 

2 Zfc Parafermions 

In this section we review the conformal properties of si{2) parafermions. In particular. 



we review the A-D-E classifications ||Tl| , |I5| of torus partition functions and present general 
expressions for the cylinder partition functions following the general framework of Behrend, 
Pearce, Petkova and Zuber ||^ on the assumption that the conformal boundary conditions 
are labelled by nodes of the tensor product graph 

G®Z2fc (2.1) 

where the Coxeter number of G is g = k + 2 and Z2fc denotes the oriented cyclic graph with 
2k nodes. 

2.1 Conformal data 

The si{2) Zfc parafermionic conformal field theories are u{l) cosets of the level k, si{2) WZW 
models 

iff i^-^' 

with central charges 

3A; 2(A;-1) , , 

c=- 1 = ^^ 2.3 

k+2 k+2 ^ ' 
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The conformal weights of primary fields (f){e,m) are given by 

+ 2) 

and the associated level k characters 



A^e,m) = /n , 0<i<k, 0<|m|<£, i-me2Z (2.4) 



X{e,m){q) = viQ) c^iq) (2.5) 



are given by the string functions [T5| 

^ 97 — r^(t m) 

(^q^ = ^ y^^ ^_lY+SqS{s+l)/2+r{r+l)/2+rs{k+l) 
^ ^'^^ r,s=0 

X |'gS(£-m)/2+r{f+m)/2 _ ^fc+l-f+s(2A:+2-£+m)/2+r(2fc+2-f-m)/2 j ^2 g) 

where 77(g) = Il^il-'- ~ 0'") Dedekind eta function and q = e^'^^'^ is the modular 

parameter. 

The conformal weights and characters are extended to the grid 

0<£<A;, 0<m<2A;-l, ^-mG 2Z (2.7) 
by identifying the underlying operators 

0(£,m) = -m) = 0(Am+2fc) = (p*k-i,k-m) (2-8) 

where * denotes conjugation. The conformal weights A(£^m) and characters X(£,m)(Q') are of 
course not effected by conjugation. The parafermionic "Kac tables" for the critical Ising 
(G = As, k = 2), critical 3-state Potts {G = A/^, k = 3) and {G = A^, k = 4) models are 
shown in Table |1|. Each operator appears twice in these grids. To remove this redundancy 
we could fix a fundamental domain to be the left half of the Kac table. But in order to stress 
the underlying structure of the fusion algebras we will sometimes make other choices, such 
as the even sub-lattice with i and m even. 

Under modular transformations, the level k parafermionic characters transform as 

Xii,m){e'n = E Sii,^)^''"''^ Xie',m'){e-'-'^') (2.9) 

{e'.m') 

where the modular matrix S satisfies 

S^ = S-\ S^ = I, S^ = C (2.10) 

and C is the matrix implementing conjugation. If {£, m) ranges over the full conformal grid 
then S = Swzw^Sj^^^ is just the tensor product of the sl{2)k Wess-Zumino-Witten modular 
matrix and the s£{2k)i modular matrix. Restricting to a fundamental domain, the entries 
of the k{k + l)/2 X k{k + l)/2 modular matrix are simply renormalized by a factor of 2 and 
given explicitly by 

c u'm') 1 • 7r(£+l)(f + l) mmm' 

SitmV ' ' = , =sm— — , ^ exp — ; (2.11) 

^ ' ^ v/A;(A; + 2) (fc + 2) k ^ ' 
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^ {k = 2) 



2 


1 

2 









1 




1 

16 




1 

16 










1 

2 





1 2 3 m e Z4 
> • > • > m > m > 



^4 {k = 3) 



3 




2 
3 









2* 
3 


2 


2 

5 




1 * 
15 




1 

15 




1 




1 

15 




2 
5 




1 * 
15 










2* 
3 




2 
3 





1 2 3 4 SmeZg 
> m > m > m > m > m > m > 



A5 (fc = 4) 



T ^ 


1 




3* 
4 









3 
4 




<• 3 




9 

16 




1 * 

16 




1 

16 




9 * 
16 


<• 2 


1 

3 








1 

3 




1 

12 




" 1 




1 




9 * 




9 




1 * 




16 




16 




16 




16 


1 







3 
4 




1 




3* 
4 





1 2 3 4 5 6 7 meZs 
> m > • > m > m > m > m > • > m > 

Table 1: Conformal grids of conformal weights for critical Ising (^43,^4), W3 {A^^Zq) and 
critical (745, Zg) models. The conjugation {l,m) 1— > {k — l,k — m) is realized as a central 
inversion symmetry through the point shown as *. 
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2.2 



parafermions on a torus 



The s£{2) parafermions on a torus have been classified by Gepner and Qiu [jTT| using 



modular invariance. They distinguish principal and non-principal cases and their results are 
summarised in Tables ^ and ^. Here we will be primarily concerned with the principal cases 
for which the identification with lattice models is clear. We conjecture that non-principal 
models are associated with fused models lETI, 



2.3 Zfc parafermions on a cylinder 

In this section we consider the parafermion models on a cylinder. We restrict our attention 
to the principal cases (G, Z2k) where G is of A, D or E type. When G is of A type, the 
theory is diagonal and the cylinder partition functions are 

Z{£i,mi)\i£2,m2){(l) = ^ ^(^,m) {^i ,mi) ^^^'""^ V(^,m) (?) (2-12) 
{e,m) 

where (£i,mi) labels the boundary condition on the left, (£2, '"^2) labels the boundary con- 
dition on the right and the fusion coefficients are given by the Verlinde formula 

Q (i',m') Q (e'.m') Q (e',m') * 

iV(£,m) (£i,mi)' = Q (Z>7) l2.idj 

(£',m') 

The parafermionic cylinder partition functions for (A4, Zg) are shown in Table ^ The char- 
acters associated to conjugate primary fields, although equal, have not been used so that 
this table encodes the complete fusion algebra with the conjugate operators shown in the 
first column. 

If G is of D or ii^ type, the boundary conditions are labelled by (a, m) G (G, Z2A;) and 
the cylinder partition functions are given by 0: 

^(ai,mi)|(a2,m2)(9) = ^(^,m) (ai ^"^'"'''''X(<?,m) (?) (2-14) 

(£,m) 

where the sum is restricted to a fundamental domain of the Kac table. The intertwiners are 

a (a'-l,m') g, ,(a',m') iTf, (a\m')* 

T/ (a2,m2) ^{e,m) ^(ai,mi) ^(a2,m2) -[^j (OTr:\ 

V(Am)(ai,mi)' =2^ g (a'-l,m') (2.15) 

(a',m') ^ ' ^ 

where the sum is over a' G Exp(G) the Coxeter exponents of G and m! G Z2A;. The eigen- 
vectors of the adjacency matrix of the graph G ® Z2A; are given by the tensor product 



(a' m'-) , a' mmm 

r' ' = i^a exp ; 

k 

where ^p"'' is the eigenvector, labelled by a' G Exp(G), of the adjacency matrix G. 
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2k-l k 



m=0 1=0 



2k-l ( fe/2-1 

)| + X(k-e,m) X{1 ,m) ,m) 



m=0 I <=() 

^cvcn 



2fc-l fc 

;, m)' 



m=0 /=o 

eeven 



2fc-i r fe fc 

(L'2„+i,Z8„_4) : IE { E IX{^,m)r+ E X(fe-^,m) 



{k = An- 2) 



m=0 I fco fco 
V £ evei 



I odd 



2fc-l 



,m) + X{6,m)r + IX(3,m) + X(7,m) ^ + IX{4,m) + X(10,m)|^} 

{k = 10) 



m=0 



2k-l 



{Et, Z32) : I ^ { |X(0,m) + X(16,m) r + IX(4,m) + X(12,m) ^ + IX(6,m) + X(10,m) ^ 

m=0 

= 16) +|X(8,m)r + (X(2,m) + X(14,m))X(8,m) + X(8,m) (X(2,m) + X(14,m))} 



2k-l 



{Es, Zse) : | ^ | |X(0,m) + X(10,m) + X(18,m) + X(28,m) r 

m=0 

{k = 28) +|X(6,m) + X(12,m) + X(16,m) + X(22,m)r} 



Table 2: The A-D-E classification of modular invariant partition functions for principal 
theories. 



k 



k 



{D2n+2,^2f3) ■ I ^ ^ {X(i,a+)+Xik-e,a+)){X{£,a-)+X(k-e,a-)) 

2a ^ even 

{k = 4n) 

{fc k 
X(£,a+)X(£,a-) + _2-^ X{k-i,a+) X{£,a- 
^ even i. odd 

(A; = 4n - 2) 

(-^6,^2/3): I E {(X(0,a+) +X(6,a+))(X(0,a-) +X(6,a-)) 

(/C = 10) + (X(3,a+) +X(7,a+))(X(3,a-) +X(7,a-)) + (X(4,a+) + X(10,a+)) (X(4,a-) + X(10,a-)) } 



{Et,Z2p) : \ E {(X(0,a+) +X(16,a+))(X(0,a-) +X(16,a-)) 

+ (X(4,a+) + X(12,«+))(X(4,a-) + X(12,a-)) + (X(6,a+) + X(10,Q+)) (X(6,Q-) + X(10,a-)) 
= 16) +X(8,a+)X(8,a-) + (X(2,a+) + X(14,a+))X(8,a-) + X(8,a+)(X(2,a-) + X(14,a 



(£^8,^2/3): i E {(X(0,a+)+X(10,a+)+X(18,a+)+X(28,a+))(X(0,a-)+X(10,a-)+X(18,a-)+X(28,a- 

xeZ2;3,2/eZi2ci 

= 28) +(X(6,a+)+X(12,a+)+X(16,a+)+X(22,a+))(X(6,a-)+X(12,a-)+X(16,a-)+X(22,a- 



Table 3: The A-D-E classification of modular invariant partition functions for non-principal 
Zfc theories. Here a 7^ 1 is a divisor of A; = afi and a'^ — ax + (3y, — ax — (5y. 
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Right 



(£,m) 


(0,0) 


(2,0) 


(0,2) 


(2,2) 


(0,4) 


(2,4) 


(0,0) 


X(o,o) 


X(0,2) 


X{0,4) 


X{2,0) 


X{2,2) 


X{2,4) 


(0,2) 


X(0,4) 


X(o,o) 


X{0,2) 


X{2,4) 


X{2,0) 


X{2,2) 


(0,4) 


X(0,2) 


X(0,4) 


X(o,o) 


X(2,2) 


X(2,4) 


X(2,0) 


(2,0) 


X(2,0) 


X(2,2) 


X(2,4) 


X(0,0) + X(2,0) 


X(0,2) + X(2,2) 


X(0,4) + X{2,4) 


(2,2) 


X(2,4) 


X(2,0) 


X(2,2) 


X{0,4) + X(2,4) 


X(0,0) + X(2,0) 


X(0,2) + X{2,2) 


(2,4) 


X(2,2) 


X(2,4), 


X(2,0) 


X(0,2) + X(2,2) 


X(0,4) + X(2,4) 


X(0,0) + X{2,0) 



Table 4: Cylinder partition functions ZL\ji for the model {A4^,Zq) ~ (T2,Z3). 

The intertwiner can alternatively be defined as the tensor product of intertwiners for 
each graph 

TA , , Aa2,m2) _ t/'^ a2 T/^ m,2 /o i yX 

>/(£,m) (ai,mi)' " l^^+l, ai Wl, mi {^■^'> 

where the intertwiners for G and Z = %2k are defined by a formulas similar to ( p.l5| ). The 
intertwiner for the graph G satisfies the recurrence relation 

y- ■■= v; ■ vt, - vt, (2.18) 



where Vi = I and = G, the adjacency matrix of G. The intertwiner for Zi2k is simply a 
power of the adjacency matrix = {1^2k)"^ of the directed graph Z2k- 

The parafermionic cylinder partition functions for (As, Zg) and (D4, Zg) are shown in 
Tables |^-^. In the case the characters associated with conjugate fields have been kept 
separated in order to encode the complete fusion algebra, but have been identified in the D4 
case to make the Z3 symmetry apparent when these (D4, Zg) partition functions are written 
in terms of the fundamental string functions. 

In the general D and E cases, these partition functions are best understood in terms 
of extended (block) characters 0. These extended characters X(a,m) are linear forms in the 
{AL,7j2k) characters involving the fundamental intertwiner 

X(a,m)(g) ■■=J2^ii,n)m^"'""^Xie,n){(l) (2-19) 

where the sum is over a fundamental domain. Note that for the (-Dfc_,_2, Z2fc) models the 
conformal dimensions of the two primary fields appearing in the extended characters differ 
by integers or half integers 



. . i{i + 2)-{k-i){k-i + 2) i-l 
A(£,,n) - ^{k~e,rn) 4:{k + 2) ~ ^ 



Note also that the parafermionic symmetries ( |2.^ ) imply that these models are in fact Z^ 
symmetric. 
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(£,m) 


(0,0) 


(0,2) 


(0,4) 


(0,6) 


(2,0) 


(2,2) 


(1,1) 


(1,3) 


(1,8) 


(1,7) 


(0,0) 




"-(0,2) 




^0,6; 


Y/o A^ 


"(2,2) 


xr^ 1 ^ 
^(1,1) 


xr^ Q^ 
^(1,3) 


X/^ K^ 
^(1,6) 


x/1 T^ 
^(1,7) 


(0,2) 






-^(0,2) 




-1(2,2) 


"(2,0) 


X/^ •7\ 
^(1,7) 


x/1 1 ^ 
^(1,1) 


xr^ Q^ 
^(1,3) 


^(1,5) 


(0,4) 




(,0,oj 


(,0,UJ 


^(0,2) 


Yz-o A^ 
^(2,0) 


"(2,2) 


Y/1 K\ 

^(1,5) 


^(1,7) 


Y/1 1 \ 
^(1,1) 


^(1,3) 




^^(0,2) 


^(0,4) 


"(0,6) 


"^(0,0) 


"(2,2) 


"(2,0) 


"(1,3) 


"(1,5) 


"(1,7) 


"(1,1) 


(2,0) 


X(2,0) 


"(2,2) 


"(2,0) 


"(2,2) 


^(0,0)+X(2,0)+X(4,0) 


>^(0,2)+>;(2,2)+X(4,2) 


^(1,1)+>^(1,5) 


'^(1,3)+'^(1,7) 


X(l,l)+X(i,5) 


X(l,3)+X(l,7) 


(2,2) 


X(2,2) 


"(2,0) 


"(2,2) 


"(2,0) 


X(0,2)+'<:(2,2)+>^(4,2) 


'<:(0,0)+X(2,0)+'^(4,0) 


X(i,3)+'<:(1,7) 


'^(1,1)+X(1,5) 


^^(1,3) +^(1,7) 


X(i,i)+'<:(1,6) 


(1,1) 


^^(1.7) 


"(1,1) 


"(1,3) 


"(1,5) 


"(1,3) +"(1,7) 


"(1,1)+"(1,5) 


"(0,0) +"(2,0) 


"(0,2)+X(2,2) 


"(0,4) +^(2,0) 


"(0,6) +"(2,2) 


(1.3) 


'^(l.S) 


^^(1,7) 


"(1,1) 


"(1,3) 


'^(i,i)+'<:(i,5) 


'<:(i,3)+>^(i,7) 


'^(0,6)+>^(2,2) 


'^(0,0)+>^(2,0) 


'^(0,2)+>^(2,2) 


'^(0,4)+X(2,0) 


(1,5) 


'^(l.S) 


"(1,5) 


"(1,7) 


"(1,1) 


X(l,3)+'^(1,7) 


'^(1,1)+'^(1,5) 


X(0,4)+X(2,0) 


'^(0,6)+'^(2,2) 


'^(0,0)+'^(2,a) 


X(0,2)+X(2,2) 


(1,7) 


'^(1,1) 


"(1,3) 


"(1,5) 


"(1,7) 


^(1,1)+'^(1,5) 


'^(1,3)+'^(1,7) 


X(0,2)+X(2,2) 


'^(0,4)+X(2,0) 


'^(0,6)+>^(2,2) 


X(0,0)+>^(2,0) 



Table 5: Cylinder partition functions Zlir for (AsjZg). The rows are labelled by the left boundary condition and the columns 
^ by the right boundary condition. The table is in fact symmetric due to conjugation symmetry. 



( a . ni ) 


(AM) 


( A -2) 


(DM) 


(B.-2) 






(-4,0) 


'^(O.O) +5^(4,0) 


2X(o,2) 


"(2,0) 


"(2,2) 


X(l,l)+X(i,3) 


X(l,l)+X(l,3) 


(-4,2) 


2X(o,2) 


x(o,o)+>i:(4,o) 


"(2,2) 


"(2,0) 


X(l,l)+X(i,3) 


X(l,l)+X(i,3) 


(B,0) 


"(2.0) 


"(2,2) 


X(0,0)+>^(4,0) 


2X(o,2) 


X(l,l)+X(l,3) 


X(l,l)+X(i,3) 


(S,2) 


"(2,2) 


"(2,0) 


2X(o,2) 


X(0,0)+X(4,0) 


X(l,l)+X(l,3) 


X(l,l)+X(i,3) 


(o,i) 


'^(1,1)+>^(1,3) 


X(l,l)+X(l,3) 


X{l,l)+X(l,3) 


X(l,l)+X(l,3) 


X(0,0)+2X(2,0)+X(4,o) 


2X(0,2)+2X(2,2) 


(0,3) 


'f(l,l)+X(l,3) 


'<:(i,i)+X(i,3) 


X(l,l)+X(l,3) 


X(l,l)+X(l,3) 


2X(0,2)+2X(2,2) 


X(0,0)+2X(2,0)+X(4,0) 



Table 6: Cylinder partition functions Zl^r for {D4, Zg). O stands for the center point of D4 and A, B for any two of the three 
other points. The table is Z3 symmetric. 



(A; = 4) 



T ^ 


1 




4 









o 
4 




" 3 




y 




1 




i 




9 * 






16 




16 




16 




16 


<• 2 






]^ 








]^ 






3 




12 




3 




12 




<• 1 




1 




9 * 




9 




1 * 






16 




16 




16 




16 


i 







3 
4 




1 




3* 
4 






7 m e Zg 



> ■ > ■ > ■ > ■ > ■ > ■ > ■ > ■ > 



D4 (A; = 4) 



1 1' 
3 ' 3 



0+1 



12 ' 12 



3_|_3> 



16 ^ 16 



1 1' 
3 ' 3 



0+1 



1 2 3 4meZ8 
> • > • > • > • > • > 



Table 7: Conformal grid for the (iP4, Zg) model resulting from the folding of the (A5, Zg) model. B, C are even, O is odd. 









X(A,0) = 


X(o,o) 


+ X(4,0) 




X(A,2) = X(0,2) + X(4,2) 








X(o,i) = 


X(l,l) + X(l,5) 




X(0,3) = X(l,7) + X(l,3) 








X(B,0) = 


X(2,0) 






X(B,2) = X(2,2) 










X(c,o) = 


X(2,0) 






X(C,2) = X(2,2) 




(a, m) 


(AO) 


(A 2) 


(5,0) 




(CO) 




(0,1) 


(0,3) 


(AO) 


X{A,0) 


X(^,2) 


X(iJ,o) 


X{B,2) 


X(c,o) 


X{C,2) 


X(o,i) 


X(0,3) 


(A 2) 


X{A,2) 


X{A,0) 


X(B,2) 


X(B,0) 


X(C,2) 


X(c,o) 


X(0,3) 


X(o,i) 


(5,0) 


X(C,0) 


X{C,2) 


X(A,0) 


X(A,2) 


X(B,0) 


X(B,2) 


X(o,i) 


X{0,3) 


(i?,2) 


X(C,2) 


X(c,o) 


X(A,2) 


X{A,0) 


X(B,2) 


X(B,0) 


X(0,3) 


X(o,i) 


(CO) 


X(B,0) 


X(S,2) 


X(c,o) 


X(C,2) 


X(A,0) 


XiA,2) 


X(o,i) 


X(0,3) 


(C^,2) 


X(B,2) 


X(s,o) 


X(C,2) 


X(c,o) 


X(A,2) 


X(A,0) 


X(0,3) 


X(o,i) 


(o,i) 


X(0,3) 


X{o,i) 


X(0,3) 


X(o,i) 


X(0,3) 


X(o,i) 


X{A,0) + X(B,0) + X(C,0) 


X(A,2) + X(B,2) + X(C,2) 


(0,3) 


X(o,i) 


X(0,3) 


X(o,i) 


X(0,3) 


X(o,i) 


X(0,3) 


X(A,2) + X(B,2) + X(C,2) 


X(A,0) + X(B,0) + X(C,0) 



Table 8: Cylinder partition functions (fusion algebra) for (D4, Zg) in terms of extended characters. A is the distinguished node. 



The set of boundary conditions (a, m) for D and E theories gives rise to an extended 
fusion algebra of conformal fields (p(a,m), namely the graph fusion algebra 

<^(a,m) X <^(a',m') = ^ A^(a,m) (a',m') V(a",m") (2-21) 
(a",m") 

with structure constants given by the Verlinde like formula 

■qj, A"-',m)xTj Aa',m)TJ, (a' ,m)* 

^ mi rT7n I ^ mo Tr7.nl ^tn.n_m.-3\ 

{a' ,m) 



where the sum on (a', m) G Exp(G') x Z2fc is over the fundamental domain. In terms of the 
extended characters, the cylinder partition functions can be written as 

Z{a,m)\ia',m'){q) = ^ N(^a,m)* {a' ,m')^"' '"^ ^ X{a" ,m")iq) (2.23) 
(a",m") 

Table ^ gives the example of the D4 case. In this form the partition functions are no longer 
explicitly Z3 symmetric. There is a distinguished node which corresponds to the identity 
operator of the fusion algebra which for D4 is {A, 0). 



3 Lattice Realizations of l^k Parafermions 

3.1 A-D-E lattice models and integrable boundaries 

The principal theories are realized as the continuum scaling limit of the critical A-D-E 
lattice models with negative spectral parameter. In these models the spin states are taken 
to be the nodes of the graph G oi A, D 01 E type. Let g be the Coxeter number of G. Then 
the bulk face weights are 



W\ 



fd 






d 




c 


















y a 


b 


a 





sin(A-M) sinwVVv^. n 

: — ^ Oac + -. — r : Om, [o.lj 

sm A sm A ipb 



where u is the spectral parameter, X = Tc/g is the crossing parameter and the weights are 
understood to vanish if the adjacency condition of A, D or E type is not satisfied along 
the four edges of the face. The crossing factors ipa are the entries of the Perron-Frobenius 
eigenvector of the adjacency matrix G. If < m < A the continuum scaling limit of these 
models describes the si{2) unitary minimal models. Otherwise, ifA — 7r/2 < u < 0, the 
continuum scaling limit describes the principal parafermions with k = g — 2. 

The A-D-E models are also integrable in the presence of a boundary [^. The integrable 
boundary conditions at a conformal point are labelled by (r, a) G {Ag_2,G)- A general 
expression for the boundary weights of the (r, a) boundary condition [ffl in terms of boundary 
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edge weights is 
d 6 



d 5 



5™ C 



h /3 




sin(^ - u) sin(^ + u + rA) sin(2M) ^ 



(3.2) 



sin(20 



sin(20 



7=1 



Here ^ is a free parameter that should be thought of as a boundary field and E™{b, c)/3-y are 
edge weights specified in Behrend and Pearce [Q. The fused adjacency matrices at level 
r which appear in the summation are simply another notation for the intertwiners = 
defined previously. These are given recursively in terms of the adjacency matrix G of the 
A-D-E graph by the si{2) fusion rules 



F' 



G, 



GF'' 



r-l 



(3.3) 



The labels a,/?, 7 are bond variables or degeneracies with P = 1, . . . S 



1 F"^ 



At u = ^ the boundary weights are independent of ^ and decompose simply in terms of 
boundary edge weights 



5™ c 



d 6 

b (3 



= X]i?^''^(6,c)^^E™(rf,c) 



57 



(3.4) 



7=1 



It is precisely at m = ^, when the triangle boundary weights reduce effectively to edge weights, 
that the integrable boundary conditions correspond @| to conformal boundary conditions. 

The A-D-E face weights and boundary weights satisfy the Yang-Baxter and boundary 
Yang-Baxter equations. This ensures integrability through commuting double row transfer 
matrices. The entries of the N column double-row transfer matrix with boundary conditions 
(ri, ai) on the left and (r2, 02) on the right are given by 



^1,6) 




di 



X—u 
\ 






u 

\ 




U 

\ 



dN-idN dpf 



Pi bo bo 



^JV-l f'JV 




(3.5) 



3.2 Finite-size corrections 

The properties of the A-D-E lattice models connect to the data of the associated conformal 
field theories through the finite-size corrections to the eigenvalues of the double row transfer 
matrices. If we write the eigenvalues of Dr-^ai\r2aJiu, (,2) as 

Dn{u) =exp{-En{u)), n = 0,1,2,... (3.6) 
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then the finite-size corrections to the energies take the form 

Stt sin "3 / c \ / 1 \ 

EM = '2Nfiu) + fr„a,\r„aM + ("^ + A„ + A;„ j + o f -J , fc„ G N (3.7) 

where f{u) is the bulk free energy, /ri,ai|r2,a2('^) is the boundary free energy, c is the central 
charge, A„ is a conformal weight and the anisotropy angle is given by 

{gu, < u < X (unitary minimal) 

(3.8) 
_2{k+i)_ A — 7r/2<ti<0 (Zfe parafermions) 

where g — k + 2 is the Coxctcr number. 

Removing the bulk and boundary contributions to the partition function on a cylinder 
leads to the conformal partition function Zj\k{q) with left and right boundaries j = (ri,ai), 
k = (r2, 02). This can be expressed as a linear form in characters 

^.•|fc(?) = E^^/^^(^) (3.9) 

i 

where i is summed over the primary fields and the integers njj*^ e N give the operator 
content. For M double rows the modular parameter is 

M 

g = exp(27rir), T = i— sini? (3.10) 

where M/N is the aspect ratio of the cylinder. 

For a system to be conformally invariant it is usually demanded that it is both isotropic 
and translationally invariant. For the unitary minimal A-D-E models with < -u < A the 
geometry is that of an isotropic square lattice when u = X/2 and = 7r/2. In this case the 
alternation of rows in the double row transfer matrix also disappears since X — u — u. So 
the conformal point occurs for 

X — u = u — ^i — ^2 — X/2 (unitary minimal models) (3-11) 

In contrast, for the parafermionic A-D-E models with A — 7i/2 < u < 0, the relevant 
geometry is not that of an isotropic square lattice. Instead, for Zjt parafermions, the relevant 
choice is M = —A with i) = 271 /k. Thus for the W^^ or hard hexagon model the geometry is 
that of the triangular lattice with -i? = 27r/3. In general it follows that A — -u = 2A so the 
alternation of rows in the double row transfer matrix persists even though the second row 
is still the transpose of the first row. Moreover, to ensure that the left and right boundary 
conditions are conformal we need to choose ^1 = 2A and ^2 — —X so the conformal point is 

X — u = ^1 = 2A, u = ^2 = —X (Zfe parafermions) (3-12) 

Notice that at this conformal point the double row transfer matrix is not left- right symmetric 
and the boundary free energies for a given boundary condition (r, a) on the left and right 
are different because they have different boundary fields ^1 = 2A and ^2 = —X. 
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3.3 Bulk free energies 

In this section we obtain the bulk free energies f{u) or equivalently the partition functions 
per site k{u) = exp(— /(«)). Two A-D-E models sharing the same Coxeter number are 
related by intertwiners so their bulk free energies are the same. Their boundary free energies 
are also related. Thus we only need to find the free energies for the models. 

For general L = k + 1, the partition function per site of the Al model satisfies the 
inversion relation 

( \ ( sin(M + A) sin(M - A) sin^^u 
.(u) .{u + A) = -^^ sinf±i(. + A) = ^'-''^ 

and the crossing symmetry 

k{u) = k{X-u). (3.14) 
The inversion relation is established by keeping just the dominant terms in the TBA and 



inversion identity hierarchies [Q, using height reversal symmetry and equating just the 
bulk terms of order 2N. The inversion relation, crossing symmetry and height reversal 
symmetry determine a unique solution which is analytic and non-zero in the analyticity 
strip Re (u) G (-^^,|), which contains the physical strip Re (u) E (— | + A, 0) for the 



regime we consider. This problem has been solved, in a different context, by Baxter p^ . 
The solution for general L can be written as an integral 

, , f~^°° . , sinhfvr — X)t sinh(7r — 3A + u)t — sinh Xt sinhfu + X)t , , 

K{u)=exp sinhwt ^ ^ ^ — ^ ^ ^- dt. (3.15) 

J_^ tsinhvrt sinh(7r — 2A)t 

The solution can also be written explicitly for L even {k odd) as 

_ sin(^ - 2A) ""^^ sin(n + {2k - 1)A) sin (^ + 2kX) 
^^""^ sin A 11 sm{u + 2kX) sin |±i(n + (2/t - 1)A) ' ^ ' 

3.4 Boundary free energies 

In this section we obtain results for the boundary free energies by extending the inversion 



relation method [^, The boundary free energy contribution comes in three different 
parts 

fr,s\r',s' = fo + fr,s + f^,s'- (3-17) 

The vacuum term /o does not depend on the boundary condition. It is free of zeros and its 

f / R 

logarithm is analytic on the physical analyticity strip. The remaining contributions fr,s for 
each boundary condition (r, s) are different on the left L and on the right R and exhibit zeros 
in the physical analyticity strip. The boundary free energies are obtained by repeating the 
analysis of the TBA and inversion identity hierarchies keeping the dominant terms as 
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in the bulk calculation but this time equating the boundary terms of order 1 in the physical 
analyticity strip Re (m) G ( '"^^^ —-^ 



2 ' 2r 

Let gf, be the RHS of the bulk inversion relation (|3.13| ). Then the vacuum functional 
equation we obtain is, with noiu) = exp(— /o(m)), 

/ \ / A \ / \ / ^ \ 4 sin^ A /o 1 o\ 

Ko{u) Ko{u + A = qb{u) qb{-- - u) ^— — — — . 3.18 

2 sm(A — 2m) sm(A — 2(u + A)) 

Hence we find the solution with the required analyticity is expressed in terms of the bulk 
free energy as 

/ \ / \ / ^ A \ 2 sin^ A /o 1 

Ko{u) = k{u) k{-- -u + X) —— — 3.19 

2 sm(A — 2u) 

For L even, one has an explicit expression 
, , sin2(M + A) 

kJu) = f 3.20 

^ ^ sin(A-2M) ^ ' 

t^' sin2(M+ (2fc + l)A) sinf±l(M + 2A;A) cos f±}(M + (2A; - 1)A) 
^ sin2(M + 2A;A) sin + (2A; - 1)A) cos^{u + 2kX) 

In this case the inversion relation admits other explicit solutions. Combining the zeros and 
poles of the two bulk free energy terms for L = 1 mod 4, we find 

sin2(u + - 1)A) sinf±i(M-A) 

sin(2M - A) sinf±i(n + - 1)A) 



«^o(m) = —77^:^ L+w,,^.L-i TTTT (3-21) 



' sin2(M + (2A;- 1)A) sm^{u + 2kX) cos + (2/c - 1)A) 

^ sin2(M + 2A;A) sin f±i(M + (2fc - 1)A) cos|^(M + 2fcA) 

In the remaining case L = 3 mod 4, the vacuum boundary free energy can only be 
expressed in terms of integrals. Let go be the RHS of the functional equation ( p.l^ ). It 
is Analytic and Non-Zero in the strip Re (m) G ( , ~-|)- Furthermore the derivative /q 
approaches a Constant when Im {u) —>■ ±oo (ANZC). Hence we can introduce the Fourier 
transforms of the derivatives 

^o(fc) ■■= /o(«)e~'"^^« (3.22) 



-oo 



^/oH = / ^o(A;)e^"dA; (3.23) 



so that ( |3.18| ) becomes 



Toim + e^') = / 'i^.-'-du (3.24) 
2nx J qo{u) 
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The solution by inverse Fourier transform gives 

= 2^/^^(^ loggo(^-^)) (3.25) 

0<Re (w))<min(A,M+-|A) 

Integrating with respect to u and taking the w integration along the vertical line Re [w) = 
e > we obtain in the limit e — > 

loggo(M) L + 1 f °° Im \og qo{u-iw) 

fo{u) = + / ■ , ... dw (3.26) 

2 TT Jo smh(L + l)w 

The boundary free energy f\r,s) arises from a boundary condition of type (r, s) with a 
boundary field It does not depend on s. An integral form analogous to ( |3.26|) can be 
derived with go being replaced by the RHS of the boundary inversion relation 

sin(M + ^) sin(M — ^) sin(M + ^ + r\) sm{u — ^ — r\) 

sin^ A 

sinf±i(»-^-(r-l)A) sinf±l(» + e + A) 



^■{r,s){u) l^{r,s){u + A) 



The solution must be analytic on the analyticity strip with the same zeros as 

sin(M + Osin(M-^-rA). (3.28) 

For r odd, this inversion relation admits the explicit solution 

sin(n + sin(n - ^ - rA) 

K(r,s){u) = r^- (3.29) 

sm A 

r-l 

-pi- sin(M + ^ + 2A;A) sin(M - ^ - (2A; - 1)A) 
^ sin(M + ^ + (2A;- 1)A) sin(M - ^ - 2A;A) 



X 



sinf±i(M + e + (2A;- 1)A) sin f±i(M - ^ - 2A;A) 



sinf±i(u + ^ + 2A;A) sin f±}(u - ^ - (2A; - 1)A) 

with zeros at m = — ^, + rX, — ,^ — vr + A, C,—n + {r + 1)A. Since these last two zeros are not 
desired, the formula is only valid when these points are outside the analyticity strip, i.e., 

"'' + ^<Re(0<^-rA. (3.30) 



2 

There is another solution when L + r is even 



sin(M + ^ — A) sin('u — ^) 
sin^ A 

sin(M + C + {2{k - 1) + r)A) sin(M - ^ + 2A;A) 



'.(r,s){u) = r^- (3.31) 

sm A 



L-r 
2 

n sin(M + ^ + (2A;-l + r)A) sin(M - ^ + (2A; - 1)A) 

sin f±i(M + ^ + (2A; - 1 + r)A) sin ^{u - ^ + {2k - 1)A) 
^ sin f±i(u + e + {2{k - 1) + r)A) sin ^{u - ^ + 2{k ~ 1)A) 
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with zeros at m = — ^, ^ + r\, ^ + (r — 1)A — vr, — ^ — (r — 2)A — vr. Clearly, this formula is 
only valid when 

- rA < Re (0 < ^ - rX. (3.32) 

In the A4 and ^5 cases, these formulas allow us to compute the boundary free energies. 
The values for ^5 are listed in Table ^ 



/o(-A/2)=logi^ 

r = 1 2 3 4 

(-A/2) = ilogf ilog3 ilogf 
(3A/2) = ilog2 log I ilogf 

Table 9: The right and left boundary free energies /^^^-^ = f(r,s){u) and Z^^^-, = /(,.^s)(A — u) of 
the A5 model at the special point u = —A/2 and A — m = 3A/2. The boundary free energies 
are independent of s. 

Some symmetries of the parafermion models are reflected in the boundary free energies. 
The left / right symmetry is broken in general but the usual height reversal symmetry r, s 1— > 
r, L + 1 — s remains. Another symmetry comes from the fact that at a conformal point 
the boundary weights factorize into edge weights. The interchange of the two independent 
sublattices is then implemented by the Kac table symmetry r, s ^— > L — r, L + 1 — s. 



4 Identifying integrable and conformal boundary con- 
ditions 

For the unitary minimal A-D-E models with < m < A there is a one-to-one corre- 
spondence [H between integrable and conformal boundary conditions labelled by (r, a) G 
{Ag^2, G) that respect the symmetry of the generalized Kac table. In contrast, for the 
parafermionic A-D-E models with X — tt/2 < u < 0, the situation is more subtle and the 
conformal labels (a, m) G (G, Z2A:) do not coincide with the natural labels (r, a) arising from 
the construction ^ of integrable boundary conditions at the conformal point. It is therefore 
necessary to make an identification of the integrable and conformal boundary conditions 
from numerical data. 

For the A type models the dominant or vacuum configuration is given by the 2k 
periodic sawtooth configuration 

a = {... ,1,2,... + ,2,1,...} (4.1) 
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It is thus clear that the vacuum (£, m) = (0, 0) on the left and the right should be identified 
with the integrable boundary condition (r, a) = (1,1) on the left and the right with the 
number of faces = mod 2k and this is confirmed by numerical computation. Next, if we 
fix the boundary condition to the vacuum {i, m) = (0, 0) on the left or right we find from 
( 2.12| ) that the cylinder partition function reduces to a single character 

Z(0,0)\{e,m)iq) = ^{£,m)|(0,0)('?) = X{£,m){q) (4.2) 

If we fix mod 2k, this allows us to identify, up to possible ambiguities of conjugation, 
the integrable boundary condition (r, s) on the left or right with the conformal labels {i, m) 
from the numerically determined cylinder partition functions. Once these identifications are 
made they can be checked for consistency against all the other cylinder partition functions 
given in ( ^.12 ). 

Proceeding in this way we find the following correspondence for the Al parafermions 
with = mod 2k 

(^-1,^-1), r + seven 
' \ {k + l-s, k-r), r + sodd. ^ ' 

This is consistent with the Kac table symmetries, (r, s) ^ {k + 1 — r, k + 2 — s) and ( p.8|) , 
of both the minimal models and parafermions. The correspondence for the ^3, A4 and 

parafermion models is listed explicitly in Table [1^. The height reversal symmetry s 
k + 2 — s, fiips the tables about their middle rows and corresponds at the level of the 
conformal algebra to a fusion with the field v'(o,fc) which is of order two. The conjugacy 
(/, m) {k — I, k — m) can be seen in these tables as a central symmetry which centre has 
been indicated for the even sublattice. The centre of the usual Kac table symmetry of the 
minimal models has been indicated too. 

For a number of faces 7^ mod 2k, the left/right symmetry of this correspondence 
is broken. Keeping the same correspondence between integrable boundary conditions (r, s) 
and fields {i, m) on the left hand side, the correspondence on the right hand side for a 
number of faces A^ = 1 mod 2k is simply given by the same table where the columns are 
shifted cyclically to the left, the first column becoming the last column after a top/bottom 
fiip as shown in Table |Tl] and so on for other mod properties. Therefore, for a number of 
faces N = k mod 2k, the correspondence table is completely fiipped about the middle row 
and we get the height reversed correspondence. This result is supported by our numerical 
computations for k = 2,3,4 summarised in Section ^. 

This correspondence can be carried over to the D and E cases. For A^ = mod k, the 
result is again given by the trivial identification of A^. with Z^, for the even sublattice and its 
fiip for the odd sublattice. The cylinder partition functions in this case are obtained from 
the A cases by intertwiners Q 

L 
s=l 
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3 
2 
1 



(0,2) 


(0,0) 


(1,1)' 


•(1,1) 


(0,0) 


(0,2) 



2 
1 




1 2 r e ^2 



(1,3) 
(2,1) 




(1,1) 
(2,3) 






(2,2) 
(1,2) 




(2,2) 
(1,2) 


(1,1) 
(2,3) 




(1,3) 
(2,1) 





> ■ > ■ > ■ > 



3 m e Z4 



^4 {k = 3) 
s £ 



4 


(0,2) 


(3,1) 


(0,0) 


1 


• 


3 




(2,4) 
(2,1) 




(1,1) 
(3,4) 




(3,1) 
(1,4) 


3 


(2,0) 


(1,1) 

• 5 


(2,2) 


1 


• 


2 


(1,3) 
(3,2) 




(3,3) 
(1,2) 




(2,3) 
(2,2) 




2 


(2,2) 


(1,1) 


(2,0) 


1 


• 


1 




(2,2) 
(2,3) 




(1,3) 
(3,2) 




(3,3) 
(1,2) 


1 


(0,0) 


(3,1) 


(0,2) 


1 


• 





(1,1) 
(3,4) 




(3,1) 
(1,4) 




(2,4) 
(2,1) 





123r 012345m 
• • • > • > • > • > • > • > • > 



A (A; = 4) 
s £ 



5 


(4,0) 


(0,2) 


(4,2) 


(0,0) 


1 


> 


4 


(1,5) 
(4,1) 




(3,5) 
(2,1) 




(1,1) 

(4,5) 




(3,1) 
(2,5) 




4 


(1,3) 


(3,1) 


(1,1) 


(3,3) 




> 


3 




(2,4) 
(3,2) 




(4,4) 
(1,2) 




(2,2) 
(3,4) 




(4,2) 
(1,4) 


3 


(2,0) 


(2,2)1 


• (2*2) 


(2,0) 




> 


2 


(1,3) 
(4,3) 




(3^3) 
(2*3) 




(1,3) 
(4,3) 




(3,3) 
(2,3) 




2 


(3,3) 


(1,1) 


(3,1) 


(1,3) 




> 


1 




(2,2) 
(3,4) 




(4,2) 
(1,4) 




(2,4) 
(3,2) 




(4,4) 
(1,2) 


1 


(0,0) 


(4,2) 


(0,2) 


(4,0) 


1 







(1,1) 
(4,5) 




(3,1) 
(2,6) 




(1,5) 
(4,1) 




(3,5) 
(2,1) 






1 


2 


3 


4 


r 









1 


2 


3 


4 


5 


6 


7 



• • • • > ■ > ■ > ■ > ■ > ■ > ■ > ■ > ■ > 



Table 10: The correspondence (r, s) 1— > {£, m) between the construction labels (r, s) e ^4^ x 
Afe+i of the integrable boundary conditions and the parafermionic conformal fields {£, m) e 
Afe+i X Z2fc and its inverse, for k — 2,3,4 and = mod 2k. The inverse is a 1:2 mapping. 
The conjugation is realized in the even sub-lattices by a central inversion symmetry through 
the point indicated by *, and the minimal model symmetry centre is indicated by •. 
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^3 (A; = 2) 

s 



{k = 3) 

s 



As (A; = 4) 

s 



o 




fo 01 




4 


(3,1) 


(0,0) 


(0,0) 


(J 


(0,2) 


(4,2) 


(0,0) 


(0,0) 


2 


(1,1) 


(1,1) 




3 


(1,1) 


(2,2) 


(2,2) 


4 


(3,1) 


(1,1) 


(3,3) 


(3,3) 


1 


(0,2) 


(0,2) 




2 


(1,1) 


(2,0) 


(2,0) 


3 


(2,2) 


(2,2) 


(2,0) 


(2,0) 




1 


2 


r 


1 


(3,1) 


(0,2) 


(0,2) 


2 


(1,1) 


(3,1) 


(1,3) 


(1,3) 












1 


2 


3 


r 1 


(1:2) 


(0.2) 


(1.0) 


(1,0) 



1 2 3 4 r 



Table 11: The correspondence (r, s) ^ {£,m) between the construction labels (r, s) of the 
integrable boundary conditions on the right and the parafermionic conformal fields {£, m) 
for the Al models with L = 3, 4, 5 and number of faces = 1 mod 2k. 

ae D4 



(a,0) 


(a,2) 


(a,2) 


(a,0) 


(0,3) 


(o,i) 


(o,i) 


(0,3) 


(A,0) 


(^,2) 


(^,2) 


(A,0) 




r e 1 



(l,a) 
(4, a) 




(3,0) 
(2,a) 






(2,0) 
(3,0) 




(l.O) 
(4,0) 


(1,^) 
(4,A) 




(3, A) 
(2, A) 





1 2 3 m e Zg 
• > • > • > • > 



Table 12: The correspondence (r, a) (a, m) between integrable boundary conditions (r, a) 
and conformal fields (a, m) for the {D4, Zg) model with the number of faces N — mod 4. 
Within these tables a — B or C. The minimal model Kac table symmetry reduces to the 
left/right fiip in the A4 x D4 table. 



B C 




O 
A 



(a,2) 


(a,2) 


(a,0) 


(a*,0) 


(o,i) 


(o,i) 


(0,3) 


(0,3) 


(^,2) 


(A,2) 


(A,0) 


(A,0) 



1 2 3 4 r e A4 
• • • • 



Table 13: The correspondence (r, a) ^ (a, m) between integrable boundary conditions (r, a) 
and conformal fields (a, m) for the (D4, Zg) model with the number of faces N = 1 mod 4. 
Note that = C and C* = S. 
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Comparing our numerical results with Table ^ allows us to obtain the correspondence given 
in Table |12| for the (D4, Zg) case with the number of faces N = mod 4: 



A4 X 

fr, a) 



X Zg 

(a, r — 1) if r + a is odd, 
(a, 4 — r) if r + a is even. 



For example, the partition function for the boundary conditions (r, a) 
left and right is: 



(4.5) 
'1, O) on the 



7D4 
^io\io 



(4.6) 



s=l 



^n|ii(?) + 2^n|i3(?) + ZtiiiM = X(o,o)(g) + 2x(2,o)(g) + X(4,o)(g). 



This result agrees with the alternative way of computing this partition function using the 
correspondence (|4.5|) and the extended fusion algebra described in Table ||: 



7D4 



il) = ^(0.3)1(0,3) (^) (4-7) 

= X{A,o) (q) + X{B,o) (q) + X{c,o) (q) = xm (q) + X(4,o) (q) + X(2,o) (q) + X(2,o) (q) ■ 



More generally, the correspondence in the D and E cases relate ( p.l9| ) and ( [4.4| ) as two 
compatible ways of computing partition functions in terms of partition functions. 

As in the case, the table for = 1 mod 4 is obtained by a shifting procedure on the 
table, namely, shift the columns of the table to the left and interchange the B and C entries 
in the last column (see Table |T^). 

Using the intertwiner definition ( [4.4|) and the correspondence formula ([4.3|) , one finds 
that this result holds for any Dk^^ and for the exceptional cases Eq, Ej and Eg: 



a, m 



(a, r — 1) if r + a is odd, 
{a,k — r) if r + a is even 



(4.8) 



where a E G, r E Ak and m G ^2^. This formula also holds for G = Ak+i provided the 
ranges of r and s are adjusted to begin at 1 rather than 0. Here the labelling of the nodes 
of the D and exceptional E graphs is as follows: 



1 2 



+2' 




(1+1)* 



1 2 


4 5 


1 2 3 


5 6 


12 3 4 


6 7 
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5 Numerical Conformal Spectra 



In this section we describe our numerical analysis which confirm the results presented in 
Section ^ As explained in the previous section, the spectra for the D and E models are 
related to the spectra of the A models by intertwiners ([4.4| ) so we need only do numerics 
on the A parafermion models. The Yang-Baxter equation, boundary Yang-Baxter equation, 
reflection symmetry, inversion relation, crossing symmetry, height reversal symmetry and 
commutation of double row transfer matrices for various number of faces and all boundary 
conditions were checked numerically. The numerics were carried out for finite size row 
transfer matrices and extrapolated to large sizes. The spectra were obtained by numerical 
diagonalization of the row transfer matrices. For the Al lattice models the maximum number 
of faces in a row that we found manageable due to machine limitations was 22, 18 and 16 
respectively for L = 3, 4, 5. 

We first studied the boundary conditions which from ([4.2|) should yield single characters, 
such as (r, 1|1,1) or (l,s|l,l) with the number of faces = mod 2k. We computed 
numerically the double row transfer matrices for these boundary conditions for increasing A^ 
and numerically diagonalized them to obtain the spectra. Allowing for the contribution of 
the bulk and boundary free energies, the central charge and the geometric factor, we fit the 
data corresponding to the largest eigenvalues of these transfer matrices, in negative integers 
powers of A^ to extract the conformal weight A. This value was compared with the entries 
in the parafermion Kac table to determine (up to conjugacy) the primary field associated 
with the given boundary condition. 

In the A^ case the level is /c = 4 and machine limitations prevent us from dealing with 
more than A^ = 16 faces. Consequently, this direct method would give us only two numbers, 
for N = 8 and A^ = 16 faces respectively. However, some symmetries can be used to obtain 
more data for extrapolation to large A^ by merging sequences for related boundary conditions 
and other mod properties as we now explain. Firstly, height reversal symmetry allows us 
to merge the sequence (r, s|l, 1) for A^ = mod 2k with the sequence (r, s|l, L) for A^ = k 
mod 2k. More generally, other sequences converge to the same character, namely (r, s\m, 1) 
for N = 2k — m mod 2k and (r, s\m,L) for N = k ~ m mod 2k. Indeed, these boundary 
conditions are all compatible with the sawtooth shaped ground state. We discovered that 
this shift and flip procedure not only applies for the vacuum but for any boundary condition 
(see Table |1T]). This method of interleaving sequences allowed us to improve dramatically 
the accuracy of our numerics. A typical extrapolated numerical estimate of a conformal 
weight agrees to 5 or 6 digit accuracy with the exact entry in the parafermion Kac table. 

In addition to producing the correct conformal weight A, the numerical spectra should, 
in the case of a single character, also reproduce the correct degeneracies of the associated 
string function X{t,m)- This means that if the eigenvalues A^^^ of the double row transfer 
matrix with A^ faces are placed in decreasing order with n = 0, 1, 2, . . . then 

Ar) = Are.p(-?^a„ + <,(l)) (5.1) 
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Energy 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Numerical 

Exact 


1.04-10-'' 




0.33331 

1 
3 


0.999502 
1 


1.33363 
4 
3 


1.33304 
4 

3 


1.99795 

2 


2.01268 

2 


2.35139 
7 
3 


2.33541 
7 
3 


2.3362 
7 
3 



Table 14: Comparison of the exact values and numerical estimates of the first ten energy lev- 
els for the parafermion model with (1, 3|1, 3) boundary conditions. The cylinder partition 
function is ^13113(9) = ^(2,0)1(2,0) (?) = X(o,o)(g) + X(2,o){q) +X(4,o)(g)- 



where a„ is an increasing integer sequence encoding the degeneracies in the character in 
terms of the modular parameter q = exp(— ^ sin-i?). The leading corrections in the o(-^) 
term are negative integer powers of so a plot of the normalized log value against 1/A^ 
becomes a straight line near zero and a polynomial fit in 1/A^ gives a good extrapolation to 
= 00. 




Figure 1: The extrapolated sequences corresponding to the first ten energy levels of the 
double row transfer matrix for the parafermion model with (1, 3|1, 3) boundary conditions. 
The horizontal axis is 1/A^ where N is the number of faces. 



After investigating all the boundary conditions leading to a single character we turned 
to the other boundary conditions and compared the numerical spectra with the predictions 
of the parafermion fusion algebra. In these cases the cylinder partition functions are sums 
of characters. As a typical example, consider the (A5, Zg) parafermion model with (1, 3|1, 3) 
boundary conditions and = mod 8. According to Table [l^, the normalized partition 
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function should converge to the sum of characters: 



^13|13(g) = ^(2,0)1(2,0) (g) = X(0,0)(g) + X(2,0)(g) + X(4,0)(g) (5.2) 

= q~^/^\l + q^'^ + q + 2q'"^ + 2q^ + 3^^/=^ + Sg^^ + o(g3)) 

so the energy levels and their degeneracies should be given by the sequence 

1 4 4 7 7 7 
{a4 = {0,-,l,-, -,2,2, 3...} (5.3) 

Table and Figure |I| show a comparison of the extrapolated sequences and the exact values. 
The accuracy in other cases is at least as good as this case. As can be seen in Figure 0, 
in extrapolating the sequences, the eigenvalues cannot simply be ordered according to their 
values at a given N . However, a given eigenvalue can be uniquely identified and tracked for 
each N by examining its pattern of zeros in the complex w-plane. 

The D and E parafermion models can be studied numerically by the same direct numer- 
ical diagonalization techniques as the A-type models. However, in these cases, it is better 
to use the intertwiner relations ( [4.4| ) to relate |2^ the spectra of the D or E model to the 
spectra of various sectors (or boundary conditions) of the associated A model with the same 
Coxeter number. In this way, for example, the spectra of the D4 model with given boundary 
conditions is related to sectors of the model. Although it is not necessary to directly 
diagonalize the double row transfer matrices of the D and E parafermion models we did 
diagonalize some cases as checks on the numerics. 



6 Discussion 

In this paper we have shown how to associate integrable lattice boundary conditions to each 
conformal boundary condition (a, m) of the A-D-E si{2) parafermion models. Moreover, 
we have shown how the conformal labels (a, m) of the string functions are related to the 
labels (r, a) that naturally arise in the construction of the integrable boundary conditions. 
In contrast to the unitary minimals models, we observe ( |3.12|) that the symmetry between 
the left and right boundaries is broken for the parafermion models due to the presence of 
an intrinsic cyclic chirality. We note also that there does not exist a gauge in which the 
local Boltzmann weights of the parafermion lattice models are all positive. Since the local 
weights should represent probabilities, this means that these models are not well defined as 
lattice models in statistical mechanics. Nevertheless, as is also the case for the non-unitary 
minimal models, these models yield well-defined conformal field theories in the continuum 
scaling limit. 

In this work we have focussed on the principal parafermion models. It would be of 
interest to extend our results to the non-principal models. We expect these models to 
correspond to fused versions [pil. El] of the lattice models studied here. 
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